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CHAPTER  1 
INTRODUCTION 


The  purpose  of  this  report  is  to  provide  a  user's 
manual  for  "The  Electromagnetic  Surface  Patch  Code"  (ESP) . 
This  code  implements  a  moment  method  solution  for  thin-wires 
and  rectangular  plates.  This  report  will  describe  the  use 
of  the  ESP  code,  rather  than  the  details  of  the  code  itself 
or  it's  theoretical  background.  These  details  can  be  found 
in: 

1.  E.  H.  Newman  and  D.  M.  Pozar,  "Electromagnetic 

modeling  of  of  composite  wire  and  surface 
geometries,"  IEEE  Trans.  Antennas  and  Propagat., 
vol .  AP-26,  pp.  784-789,  Nov.  1978. 

2.  E.  H.  Newman  and  D.  M.  Pozar,  "Considerations 

for  efficient  wire/surface  modeling,"  IEEE  Trans. 
Antennas  and  Propagat.,  vol.  AP-28,  pp.  121-125, 
Jan.  1980. 


The  code  is  capable  of  treating: 

1.  thin  wires 

2.  rectagular  plates 

3.  wire  to  plate  junctions  atleast  0.1  wavelengths 
from  an  edge 

4.  plate  to  plate  junctions,  including  several  plates 
intersecting  at  a  common  edge;  the  sides  of 
intersecting  plates  need  not  "line  up" 


5.  excitation  by  delta-gap  generator  or  plane  wave. 
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6.  open  or  closed  surfaces 

The  code  computes:  j 

I 

1.  currents 

2.  input  impedance  and  admittance 

3.  efficiency  ^ 

4.  far-zone  patterns  (both  polarizations) 

I,  ] 

5.  back  or  bistatic  scattering  (theta,  phi  and  cross 
polarization) 


Some  special  features  of  the  code  are: 


1. 

simplified  input  of  problem 

computation  desired 

geometry 

and 

2. 

overlap  rodes  enforce  continuity 
plate  to  plate  junctions 

of 

current 

at 

3. 

attachment  modes  enforce  continuity 
wire  to  plate  junctions 

of 

current 

at 

4. 

a  wire  may  contact  a  plate  anywhere 
i.e.,  the  junction  need  not  have 

in 

any 

it's  area, 
relation  to 

the  plate  modes 

5.  for  a  series  of  problems  where  the  geometry  changes 
only  slightly  from  one  run  to  the  next,  only  that 
portion  of  the  impedance  matrix  which  has  changed 
need  be  computed  from  one  run  to  the  next 
(resulting  in  a  savings  in  time) 

6.  plots  of  patterns  and  wire/plate  geometry. 


Chapter  2  describes  the  program  input.  Chapter  3 
describes  the  program  output  and  contains  several  design 
examples.  Chapter  3  describes  the  various  array  dimensions, 
and  what  to  do  if  dimensions  must  be  changed.  Chapter  3 
also  describes  various  files  used  by  the  code.  A 
description  of  the  magnetic  tape,  on  which  the  code  is 
supplied,  is  given. 
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This  manual  and  the  code  it  describes  are  not 
necessarily  in  final  form.  For  example  the  author  has 
obtained  solutions  to  the  problem  of  wires  attached  near  the 
edge  of  a  plate  and  to  treating  nonrectangular  plates.  It 
is  anticipated,  that  in  the  future,  these  features  will  be 
integrated  into  the  code.  Another  possibilitie  for  a  new 
feature  includes  a  flat  lossy  earth.  Also,  as  feedback  is 
obtained  from  users,  the  code  may  change  to  reflect  their 
suggestions. 

A  considerable  effort  has  been  made  to  validate  the 
accuracy  of  the  ESP  code.  Computed  results  have  been 
checked  against  measurements  and  other  solutions  (wire-grid 
solutions  and  high  frequency  solutions  such  the  geometrical 
theory  of  diffraction).  However,  experience  has  shown  that 
with  complex  codes  there  is  always  a  combination  of  inputs 
or  a  geometry  for  which  the  code  will  fail.  When  writing 
and  testing  the  code  we  generally  had  in  mind  the  wires 
representing  antennas  and  the  plates  representing  support 
structure  or  other  scattering  obstacles.  We  have  not  tested 
such  things  as  wires  forming  transmission  lines,  plates 
forming  waveguides  (we  have  successfully  modeled  a  cavity 
antenna  with  plates),  or  two  closely  spaced  plates  forming  a 
thick  plate.  We  have  modeled  microstrip  antennas  with 
plates,  but  this  required  a  special  purpose  code.  The  point 
is  that  the  results  of  this  or  any  other  code  should  not  be 
trusted  absolutely.  Atleast  some  of  the  results  should  be 
compared  with  measurements  or  other  computations  on  a 
similar  problem.  This  is  especially  true  if  the  problem  is 
considerably  different  from  those  previously  run. 
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PROGRAM  INPUT 


2.1  INTRODUCTION 

This  chapter  describes  the  input  to  the  Electromagnetic 
Surface  Patch  Code  (ESP).  That  is,  this  chapter  will  give 
the  method  for  describing  to  the  program  the  problem 
geometry  and  the  type  of  computation  desired.  This  input 
data  can  be  broken  into  four  parts  as  follows: 

1.  the  wire  geometry 

2.  the  plate  geometry 

3.  the  type  of  computation  desired 

4.  miscellanous  parameters. 

Two  methods  of  inputting  this  information  are  employed.  The 
wire  geometry  and  a  description  of  wire/plate  junctions  are 
input  by  the  user  writing  a  FORTRAN  subroutine  called  WGEOM, 
or  by  an  input  file  (logical  unit  5) .  All  other  quantities 
are  read  in  via  the  input  file  (logical  unit  5) .  To  code 
user's  who  are  used  to  inputting  all  data  via  an  input  file, 
the  use  of  a  FORTRAN  subroutine  to  describe  the  wire 
geometry  may  seem  like  an  unnecessary  complication. 
However,  we  have  found  that  it  often  results  in  an  easier, 
more  versatile,  faster,  and  more  reliable  method  of 
inputting  the  wire  geometry.  The  reason  is  that  the  wires 
usually  represent  antennas  or  other  structures  which  have 
geometrically  a  very  regular,  periodic,  or  simple  shape, 
expressable  in  mathematical  terms.  Good  examples  are  linear 
monopoles  or  dipoles,  helix  or  loop  antennas,  log  periodic 
antennas,  and  antenna  arrays.  In  any  case,  the  user  has  a 
choice  of  an  input  file  or  the  subroutine. 
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2.2  INPUT  READ  STATEMENTS 

The  details  of  the  input  data  will  now  be  presented. 
Figure  1  shows  the  15  input  or  READ  statements.  The  * 
implies  free  format  read.  DO  LOOP'S  are  shown  which 
indicate  the  order  of  execution. 


2.2.1  READ  1 

READ  1  inputs  several  run  control  parameters  as 

follows: 

NGO  =  0  implies  input  and  printout  geometry  and  then  stop 
(i.e.,  do  not  perform  computations). 

=  1  implies  execute  entire  run. 

NPRINT  =  0  implies  printout  certain  input  parameters 

(frequency,  wire  radius,  integration  parameters,  etc.). 

=  1  implies  print  wire  and  plate  geometry. 

=  2  implies  print  both. 

=  3  implies  print  neither. 

NRUNS  =  the  number  of  runs,  i.e.,  the  limit  of  the  DO  700  loop 
in  Figure  1. 

NWGS  =  for  each  run,  the  number  of  wire  geometries,  i.e., 
the  limit  of  the  DO  600  loop  in  Figure  1. 

IWR  =  0  implies  do  not  print  solution  or  current  vector. 

=  1  implies  print  solution  or  current  vector,  as  well 
as  detailed  geometry  for  wire  and  plate  mode  locations. 

IWRZT  =  0  implies  do  not  print  impedance  matrix. 

=  1  implies  print  impedance  matrix. 

INT  =  the  number  of  Simpson  rule  intervals  for  numerical 
integrations  on  wire  modes.  Typically  chosen  as  4  and 
always  an  even  number.  If  INT  =  0  expressions 
involving  exponential  integrals  are  used  for  wire  to  wire 
to  wire  impedances.  These  closed  form  expressions 
require  more  computer  time  than  a  numerical  integration 
with  INT  =  4.  Self  or  overlapping  impedances  are  always 
done  with  the  closed  form  expressions,  which  are  more 
accurate  than  numerical  integration. 

INTP  =  number  of  intervals  for  numerical  integration  on 
plate  modes.  Typically  chosen  as  10,  and  always  an 
even  number. 

INTD  =  the  number  of  intervals  for  numerical  integration  on 
disk  component  of  attachment  modes.  Typically  chosen  as 
18,  and  always  an  even  number. 

INWR  =  0  implies  geometry  does  not  contain  any  wires. 

=  1  implies  geometry  does  contain  wires. 

IRGM  =  0  implies  that  the  wire  and  attachment  geometry  is 
to  be  generated  in  subroutine  WGEOM. 

=  1  implies  that  the  wire  and  attachment  geometry  is  to  be 
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READ  (  5  ,  *)  IXIO.NPRINT,  NRUNS,NWGS,  IWR,  IWRZT,  INT,  IMTP,  INTD,  INWR,  IRGM  (1) 
READ  (  5  ,  *)  IFE,  I  PEE,  NDFE,  PHPE  (2) 

READ  (5  ,  ♦)  I  FA,  I  PEA,  NDFA.THFA  (3) 

READ  (  5  ,  *  )  ISE,  IPSL.NDSE,  PHSE , TH  IN ,  PH  IN  (  4) 

READ (  5  ,  *) ISA, IPSA, NDS  A , TH  SA  (5) 

D0700NRUN=1 ,NRUNS 

READ ( 5 , * ) FMC , CMM , A  (6) 

RE  AD  (5,*)  NPLTS,  IOVL  (7) 

I F ( NPLTS . EQ. 0 ) G0T046  0 
DO420MPL=l , NPLTS 

READ  ( 5  ,  *  )  Kill  2N(  NPL )  ,NM23N(MPL)  ,IPN(MPI,)  (8) 

D04  20  NCNR-l ,3 

READ  ( 5  ,  * )  PCD  ( 1 ,  NCNR,  NPL)  ,  PCN(  2  ,  NCNR,  NPL)  ,  PCN  (  3  ,  NCHR,  NPL)  (9) 

420  CONTINUE 
460  CONTINUE 

DO  600  NWG=1,NWGS 

READ  ( 5  ,  * )  IWRZM ,  IRDZM  (10) 

IF(IMWR.EO.O) GCT02773 
IF (IRGM. EQ.0) GGT02800 

READ  (  5  ,  *  )  NM,  MI>,  NAT,  NFPT  (11) 

D02  8101  =  1  , NP 

READ(5,*)X(I)  , Y  ( I )  ,Z(I)  (12) 

2810  CONTINUE 

D02  82  0  1  =  1  ,  NM 

READ(5,*)IA(I) ,IB(I)  (13) 

2820  CONTINUE 

DO2830 1=1 , NFPT 

IF  ( NFPT.  GE.  1 )  READ  ( 5  ,  *)  I FH ,  IAB,VLG,ZL  (14) 

2830  CONTINUE 

IF ( NAT . EQ. 0 ) GOT02  850 
D02  840 1  =  1  , NAT 

READ (5,*)  NAS,  IAB,  NPLA(I )  ,VCA(I)  , ZLDA  (I)  ,BDSK(I  )  (15) 


2840  CONTINUE 
GOTO2850 

2800  CALLWGEOM  (IA,  IB,  X,  Y,  Z,  NM,  NP,  NAT,  MSA,  NPLA,  VGA,  BDSK, 
2  ZLDA ,  NWG  ,  VG,  ZLD,l-V) 

2850  CONTINUE 
2773  CONTINUE 


***  MAIN  BODY  OF  PROGRAM 


600  CONTINur. 
700  CONTINUE 


Figure  1  -  The  15  input  or  READ  statements 
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read  in  via  the  input  file. 


2.2.2  READ'S  2-5 

READ'S  2-5  specify  the  far-zone  patterns  desired. 
READ'S  2  and  3  are  for  elevation  and  azimuth  radiation 
patterns,  respectivily ,  while  4  and  5  are  for  elevation  and 
azimuth  scattering,  respectively.  Specifically: 

IFE  =  0  implies  do  not  compute  far-zone  radiation  pattern 
in  the  elevation  plane. 

=  1  implies  compute  pattern. 

IPFE  =  0  implies  do  not  plot  far-zone  radiation  pattern  in 
the  elevation  plane. 

=  1  implies  plot  pattern. 

NDFE  =  angle  increment  in  degrees  for  far-zone  radiation 
pattern  in  the  elevation  plane  (should  be  evenly 
divisible  into  360) . 

PHFE  =  phi  angle  in  degrees  for  far  zone  radiation  pattern 
in  the  elevation  plane. 

I FA,  IPFA,  and  NDFA  =  same  as  IFE,  IPFE,  and  NDFE  for  far- 
zone  radiation  pattern  in  the  azimuth  plane. 

THFA  =  theta  angle  in  degrees  for  far-zone  radiation  pattern 
in  the  azimuth  plane. 

ISE  =  0  implies  do  not  compute  far-zone  scattering  pattern 
in  the  elevation  plane 

=  1  implies  compute  backscatter  pattern 
=  2  implies  compute  bistatic  scatter  pattern. 

IPSE,  NDSE,  and  PHSE  =  same  as  I*'FE,  NDFE,  and  PHFE  except 
that  they  are  for  far  zone  scattering. 

THIN  and  PHIN  =  theta  and  phi  direction  of  the  incident  wave 
for  all  bistatic  scattering  (i.e.,  ISE  or  ISA  =  2). 

ISA,  IPSA,  and  NDSA  =  same  as  ISE,  IPSE,  and  NDSE  except 
that  they  are  for  scattering  in  the  azimuth  plane. 

THSA  =  theta  angle  in  degrees  for  scattering  pattern  in 
the  azimuth  plane. 

If  ISE  or  ISA  are  set  to  -1  or  -2,  then  it  will  have 
the  same  effect  as  setting  them  to  1  or  2,  respectively, 
except  that  an  incident  image  wave  will  be  included.  That 
is,  if  a  theta  polarized  wave  is  incident  from  (theta, phi) , 
then  a  theta  polarized  wave  from  (pi-theta, phi)  will  be 
included.  If  a  phi  polarized  wave  is  incident  from 
(theta, phi) ,  then  a  -phi  polarized  wave  from  (pi-theta, phi) 
will  be  included.  The  image  plane  is  the  xy  plane.  This 

option  is  of  use  in  treating  problems  over  an  infinite 

ground  plane  using  image  theory.  Note  that  this  option 

automatically  inserts  the  image  wave,  however,  it  is  the 

user's  responsibility  to  insert  the  image  of  the  wire/plate 
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geometry . 

Note  that  one  can  not  mix  patterns  on  the  same  run. 
Thus,  on  the  same  run,  one  can  not  get  a  radiation  pattern 
and  a  scattering  pattern,  or  a  back  and  bistatic  scattering 
pattern.  Also  it  is  not  possible,  on  the  same  run,  to  get 
two  or  more  of  the  same  pattern  type.  Thus,  one  may  obtain 
an  azimuth  and  a  elevation  plane  pattern  on  the  same  run, 
but  not  two  azimuth  or  two  elevation  plane  patterns.  For 
the  user  who  wishes  a  combination  of  patterns  not  permitted 
on  a  single  run,  READ  10  will  permit  any  combination  of 
patterns  to  be  effeciently  obtained  on  several  successive 
runs . 


2.2.3  READ  6 

READ  6  inputs  the  following: 

FMC  =  frequency  in  MHZ. 

CMM  =  wire  conductivity  in  megamho/meter. 

=  -1.0  implies  infinite  conductivity. 
A  =  wire  radius  in  meters. 


2.2.4  READ'S  7-9 

READ'S  7-9  input  the  plate  geometry.  READ  7  inputs: 

NPLTS  =  the  number  of  rectangular  plates. 

IOVL  =  indicator  for  overlap  modes. 

=  0  implies  no  overlap  modes. 

=  1  implies  insert  overlap  modes  only  if  one  (or  both) 
of  the  contacting  plates  has  a  current  polarization 
mode  perpendicular  to  the  common  edge. 

=  2  implies  insert  overlap  modes  whenever  two  or  more 
plates  intersect  along  a  common  edge.  The  corners  of  the 
intersecting  plates  need  not  coincide.  See  Design  Example 
V  in  the  next  section. 

One  feature  of  the  code  is  that  it  automatically  checks 
for  plates  which  intersect  along  a  common  edge.  Surface 
patch  dipole,  or  overlap,  or  hinge  modes  can  then  be 
inserted  to  enforce  continuity  of  the  transverse  current 
across  the  common  edge.  Inserting  these  overlap  modes  is 
equivalent  to  saying  that  the  two  (or  more)  plates  are  in 
physical  contact.  Not  inserting  the  overlap  modes  is  saying 
that  the  two  (or  more)  plates  are  close  but  do  not  touch. 
The  parameter  IOVL  permits  the  user  to  select  which  of  these 
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cases  is  appropriate  to  his  problem.  Plates  are  considered 
to  have  a  common  edge  if  the  edges  are  within  TOUCH  of 
coinciding.  TOUCH  is  set  in  the  main  program  at  0.001 
wavelengths.  When  M  plates  are  detected  as  having  a  common 
edge,  M-l  sets  of  overlap  modes  are  inserted  (assuming  IOVL 
was  properly  set) .  For  most  applications  the  user  will  set 
IOVL  =  2. 

For  each  of  the  NPLTS  plates,  READ  8  is  executed  once 
and  READ  9  is  executed  three  times.  Together  these  two 
READ'S  specify  the  location,  segmentation,  and  mode  current 
polarization  on  a  given  plate.  If  NPLTS  =  0,  then  the  code 
skips  READ'S  8  and  9. 

Figure  2  shows  a  typical  rectangular  plate.  The 
location  of  the  plate  is  defined  by  the  coordinates  of  any 
three  consecutive  (clockwise  or  counterclockwise)  corners  of 
the  plate.  For  the  purpose  of  placing  surface  patch  dipole 
modes  on  the  plate,  the  plate  must  be  divided  into  smaller 
rectangular  segments  or  monopoles.  The  user  specifies  the 
number  of  segments  from  point  1  to  point  2  (=  5  in  Figure  2) 
and  the  number  of  segments  from  point  2  to  point  3  (=  4  in 
Figure  2) .  In  choosing  segment  size  the  general  rule  is 
that  the  segment  length  or  width  should  not  exceed  a  quarter 
wavelength.  Our  experience  indicates  that  a  segment  size  of 
a  quarter  wavelength  yields  reasonable  results  for  most 
applications.  If  more  accuracy  is  required,  one  can  try 
reducing  the  segment  size  to  0.2  or  0.15  wavelengths. 
However,  this  is  as  small  a  segmentation  that  one  would 
ordinarily  use.  If  one  uses  the  quarter  wavelength 
segments,  then  there  will  be  12  surface  patch  dipole  modes 
per  square  wavelength  of  surface  per  current  polarization, 
or,  24  total  modes  per  square  wavelength  of  surface.  If  the 
user  specifies  NM12  segments  in  the  1  to  2  direction  and 
NM23  segments  in  the  2  to  3  direction,  then  (assuming  the 
user  has  specified  two  orthogonal  current  polarizations  on 
this  plate)  there  will  be  NM23 * (NM12-1)  modes  polarized  in 
the  1  to  2  direction,  and  NM12* (NM23-1)  modes  polarized  in 
the  2  to  3  direction.  The  plate  in  Figure  2  would  have  16 
modes  in  the  1  to  2  direction  and  15  modes  in  the  2  to  3 
direction,  for  a  total  of  31  modes.  Specifically,  READ  8 
inputs : 

NM12(NPL)  =  the  number  of  segments  on  the  NPL  plate  in 
the  1  to  2  direction. 

NM23 (NPL)  =  the  number  of  segments  on  the  NPL  plate  in 
the  2  to  3  direction. 

IPN(NPL)  =  1  implies  place  modes  on  the  NPL  plate  with 
polarization  in  the  1  to  2  direction  only. 

*  2  implies  place  modes  on  the  NPL  plate  with 
polarization  in  the  2  to  3  direction  only. 
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=  3  implies  place  modes  on  the  NPL  plate  with 
polarization  in  the  1  to  2  and  2  to  3  direction. 

=  0  implies  place  no  modes  on  the  NPL  plate.  If  this  option 
is  chosen,  then  the  only  way  that  current  can  flow  on  the 
NPL  plate  is  if  it  contacts  another  plate  and  overlap  modes 
are  placed  on  it. 

The  above  method  for  specifying  the  modal  layout  on  a 
plate  is  ideal  in  that  it  gives  the  user  virtually  complete 
control  of,  and  yet  divorces  him  from  the  details  of  the 
overlapping  dipole  mode  layout.  Thus,  even  if  a  plate  has 
hundreds  of  modes,  the  user  need  only  specify  the  three 
integers  in  READ  8. 

Usually  one  sets  IPN(NPL)  =3  to  obtain  both 
polarizations.  However,  there  are  cases  where  the  user  may 
know  that  one  of  the  polarizations  is  not  significant.  In 
this  case,  one  can  save  computer  time  and  storage  by  setting 
IPN(NPL)  =  1  or  2. 

READ  9  is  executed  three  times  for  each  plate,  each 
time  inputting  three  floatting  point  numbers.  The  first 
time  one  inputs  the  x,  y,  and  z  coordinates  in  meters  of 
point  1  of  the  plate.  Next  one  inputs  the  coordinates  of 
point  2,  and  finally  one  inputs  the  coordinates  of  point  3. 
Specifically  READ  9  inputs: 

PCN ( I , NCNR, NPL)  =  the  I-th  coordinate  (I  =  1,2,3  for  x,y,z 
respectively)  of  the  NCNR  corner  (NCNR  =  1,2,3)  of  the 
NPL  (NPL  =  1,2,. . .NPLTS)  plate. 

As  noted  above,  READ  9  is  executed  three  times  for  each 
plate.  Thus,  for  each  plate,  one  inputs  the  three  lines: 

the  x,y,z  coordinates  of  corner  1 
the  x,y,z  coordinates  of  corner  2 
the  x,y,z  coordinates  of  corner  3. 


2.2.5  READ  10 

At  times  a  user  may  wish  to  run  several  problems  in 
succession  for  which  the  impedance  matrix  does  not  change  at 
all  or  changes  only  slightly.  For  example,  the  impedance 
matrix  will  not  change  at  all  if  one: 

1.  changes  the  far-zone  pattern  desired 


2-6 


'W!  .*■*>  f  .1  ■*  •* 


PROGRAM  INPUT 


2.  changes  the  voltage  excitation 

3.  changes  the  angle  of  incidence  in  a  bistatic 
scattering  calculation. 

Obviosly  in  these  cases  it  would  be  extremely  wasteful  to 
recompute  the  impedance  matrix.  At  other  times  the  geometry 
may  change  only  slightly  from  one  run  to  the  next.  For 
example,  suppose  one  is  faced  with  the  problem  of  locating  a 
monopole  on  a  ship  in  order  to  achieve  a  desired  impedance 
and/or  pattern.  In  order  to  solve  this  problem  one  would 
construct  a  model  of  the  ship  from  several  intersecting 
plates,  possibly  requiring  hundreds  of  surface  patch  modes. 
The  monopole  would  be  modeled  by  one  or  two  wire  modes,  and 
one  attachment  mode  would  be  required  where  the  monopole 
contacted  a  plate.  The  user  would  then  analyze  this 
configuration  for  many  monopole  locations,  looking  for  the 
one  which  best  met  his  design  goals.  The  moment  method 
impedance  matrix  for  this  (or  in  general  any)  problem  can  be 
symbolically  shown  as 


[Z] 


W/W 


W/P 


W/A 


P/W 


P/P 


P/A 


A/W 


A/P 


A/A 


where : 

W  =  wire  mode 
P  =  plate  mode 
A  =  attachment  mode. 

In  showing  IZ) ,  we  have  indicated  that  there  are,  in 
this  example,  many  more  plate  modes  than  wire  and  attachment 
modes.  Note  that  only  the  lower  triangular  part  of  the 
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symmetric  impedance  matrix  is  computed  and  stored,  the 
point  is  that  as  the  monopole  is  moved  the  P/P  block  [Z)  is 
unchanged,  since  the  number  and  location  of  the  plate  modes 
do  not  change.  Only  the  blocks  involving  wire  or  attachment 
modes  change.  Thus  a  very  considerable  time  saving  will 
result  if  on  the  first  run  the  entire  matrix  is  computed  and 
stored  (on  a  disk) .  On  subsequent  runs  the  stored  matrix  is 
read  in,  and  only  those  blocks  involving  wire  or  attachment 
modes  (i.e.,  the  lower  triangular  part  except  for  the  P/P 
block)  are  recomputed.  It  is  important  to  note  that  on  a 
problem  such  as  the  monopole  on  a  ship  considered  here  the 
P/P  block  would  involve  the  majority  of  the  computer  time. 
By  specifying  IWRZM  and  IRDZM  one  can  easily  obtain  this 
savings. 

Specifically  READ  10  inputs: 

IWRZM  =  indicator  to  write  impedance  matrix  on  a  disk 
file,  logical  unit  1. 

=  0  implies  do  not  write  out  impedance  matrix. 

=  1  implies  write  out  impedance  matrix. 

IRDZM  =  indicator  to  read  impedance  matrix  from  disk  file. 

=  0  implies  do  not  read  matrix  and  compute  entire 
new  matrix. 

=  1  implies  read  in  matrix  and  compute  new  matrix 
except  for  the  W/W  and  A/A  block. 

=  2  implies  read  in  matrix  and  compute  new  matrix 
except  for  the  P/P  block. 

=  3  implies  read  in  matrix  and  use  as  new  matrix, 
i.e.,  do  not  compute  entire  new  matrix. 

Thus,  one  would  set  IRDZM  =  2  if  the  plate  geometry  is 
unchanged  from  the  last  run,  IRDZM  =  1  .if  the  wire  and 
attachment  geometry  is  unchanged  from  the  last  run,  and 
IRDZM  =  3  if  the  entire  geometry  is  unchanged.  If  IRDZM  is 
set  >  0  then  it  is  essential  that: 

1.  IWRZM  =  1  on  a  previous  run  and 

2.  the  number  of  wire  modes,  the  number  of  plate 
modes,  and  the  number  of  attachment  modes  is 
unchanged  from  the  the  run  where  IWRZM  =  1. 


The  impedance  matrix  is  read  from  or  written  to  the  disk 
file  Z MAT. DAT. 1  on  logical  unit  1.  It  is  an  unformatted 
READ  or  WRITE. 
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2.2.6  READ'S  11-15 

READ'S  11-15  input  the  wire  geometry,  including  loads, 
generators,  and  wire  to  plate  attachments.  READ'S  11-15  are 
executed  only  if  INWR  =  1  (see  READ  1)  and  IRGM  =  1  (see 

READ  1) .  The  wire  geometry  consists  of  a  series  of 
interconnected  straight  wire  segments  or  monopoles.  Segment 
lengths  should  not  exceed  a  quarter  wavelength,  and  no  two 
intersecting  segments  should  form  an  acute  angle  less  than 
about  30  degrees.  Also,  a  single  isolated  segment  is  not 
permitted.  The  wire  geometry  input  will  be  described  with 
the  aid  of  the  example  shown  in  Figure  3.  The  structure 
consists  of  a  T  shaped  wire  with  one  load  and  one  generator. 
The  wire  consists  of  a  number  of  points,  shown  as  heavy 
black  dots  in  Figure  3,  and  segments.  An  arbitrary 
numbering  scheme  is  established  for  the  points  and  segments. 
In  Figure  3  the  point  numbers  are  shown  adjacent  to  the 
dots,  and  the  segment  numbers  are  shown  circled  next  to  the 
segments . 

READ  11  inputs  the  following: 

NM  =  the  total  number  of  segments  on  the  wire  structure. 

NP  =  the  total  number  of  points  on  the  wire  structure. 

NAT  =  the  total  number  of  wire  to  plate  attachment  points. 
NFPT  =  the  total  number  of  feed  locations  in  the  wire, 
specifically  excluding  feeds  at  wire  to  plate 
attachments.  Here  a  feed  location  is  a  location  at 
which  there  is  either  a  generator  or  a  load. 

For  the  geometry  of  Figure  3,  the  input  for  READ  11  would 
be: 

3  4  11  . 

READ  12  requires  NP  lines  of  input,  specifying  on  the 
I-th  line  the  x,y,z  coordinates  in  meters  of  the  1-th  point. 
Specifically  READ  12  inputs: 

X ( I )  =  the  x  coordinate  of  point  I  in  meters. 

Y ( I )  =  the  y  coordinate  of  point  I  in  meters. 

Z ( I )  =  the  z  coordinate  of  point  I  in  meters. 

For  the  geometry  of  Figure  3  the  input  for  READ  12  would  be 

the  NP  =  4  lines: 

0.0  0.0  0.0 
0.0  0.0  0.25 
0.0  0.0  0.5 
-0.3  0.0  0.25  . 
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READ  13  requires  NM  input  lines,  specifying  the 
endpoints  of  the  NM  segments.  Each  segment  has  two  end 
points  denoted  A  and  B.  It  is  arbitrary  which  end  the  user 
selects  as  A  and  which  as  B.  READ  13  inputs: 

IA ( I)  =  point  number  of  the  A  end  of  segment  I. 

IB (I)  =  point  number  of  the  B  end  of  segment  I. 

Arbitrarily  choosing  the  smaller  point  as  endpoint  A,  the  NM 
=  3  lines  of  input  for  READ  13  would  be: 

1  2 
2  3 

2  4 

Note  that  there  is  no  limit  to  the  number  of  wires  which  can 
intersect  at  a  given  point  (except  the  acute  angle 
limitation  described  above) . 

READ  14  inputs,  for  each  of  the  NFPT  feed  locations  in 
the  wire,  the  feed  location  and  the  complex  value  of  the 
generator  and  load  at  that  location.  In  this  code  we  always 
think  of  generators  and  loads  as  being  inserted  into 
segments,  either  by  point  A  or  by  point  B  of  the  segment. 
We  do  not  think  of  feeds  as  being  at  a  point  in  the  wire. 
For  example,  for  the  geometry  of  Figure  3,  it  is  not 
sufficient  to  say  that  the  50  ohm  load  is  "at  point  2". 
There  are  three  locations  (although  physically  close, 
electrically  very  different)  which  could  be  taken  as  "at 
point  2".  They  are  by  point  B  of  segment  1,  by  point  A  of 
segment  2,  or  (in  this  case  the  correct  choice)  by  point  A 
of  segment  3.  READ  14  inputs  for  each  of  the  NFPT  feed 
locations : 

IFM  =  segment  number  of  feed  location. 

IAB  =  0  implies  feed  is  by  point  A  of  segment  IFM. 

=  1  implies  feed  is  by  point  B  of  segment  IFM. 

VLG  =  complex  voltage  of  generator  (volts).  Note  the 
polarity  is  from  point  A  to  point  B. 

ZL  =  complex  impedance  of  the  load  (ohms) . 

For  the  geometry  of  Figure  3,  the  NFPT  =  1  line  of  input  for 
READ  14  would  be: 

3  0  (0.0, 0.0)  (50.0,0.0) 

Note  that  since  there  is  no  voltage  generator  at  the  wire 
feed  location  we  specify  zero  for  VLG. 


2-10 


PROGRAM  INPUT 


READ  15  specifies  the  wire  to  plate  attachments  and  the 
complex  values  of  the  generators  and  loads  at  the  attachment 
locations.  Specifically,  for  each  of  the  NAT  attachments  we 
input : 

NAS  =  the  segment  number  of  the  segment  which  contacts 
or  attachs  to  a  plate. 

IAB  =  0  implies  the  attachment  point  is  by  endpoint  A 
of  segment  MAS. 

=  1  implies  the  attachment  point  is  by  endpoint  B  of 
segment  NAS. 

NPLA(I)  =  plate  number  for  the  I-th  attachment. 

VGA ( I )  =  complex  generator  voltage  (volts)  at  I-th 
attachment.  Note  that  the  polarity  is  minus  at 
the  attachment  point. 

ZLDA(I)  =  complex  load  impedance  (ohms)  at  the  I-th 
attachment  point. 

BDSK(I)  =  outer  disk  radius  in  meters  to  be  used  for 
the  I-th  attachment  mode.  The  disk  is  a  circle  in 
the  plane  of  the  NPLA(I)  plate  and  centered  at  the 
I-th  attachment  point.  Two  considerations  dictate 
the  choice  of  the  disk  radius.  First,  experience 
has  shown  that  reasonably  accurate  impedance  and 
patterns  can  be  obtained  with  disk  radii  between  0.1 
and  0.25  wavelengths.  A  good  average  choice  is  0.2 
wavelengths.  Second,  the  disk  must  be  entirely 
within  the  boundaries  of  the  NPLA(I)  plate.  Thus, 
if  the  minimum  distance  from  the  attachment  point 
to  a  side  of  the  plate  is  d,  then  BDSK(I)  must  be 
chosen  less  than  d.  If  d  <  0.1  wavelengths,  then 
we  must  violate  the  first  condition.  Choosing 
BDSK(I)  <  0.1  wavelengths  will  result  in  inaccurate 
input  impedance,  but  probably  reasonably  accurate 
far-zone  patterns.  Thus  we  say  that  attachment 
points  must  be  atleast  0.1  wavelengths  from  an  edge 
if  accurate  input  impedance  data  is  required. 

Assuming  a  frequency  of  300  MHZ  (1  meter  wavelength) ,  READ 
15  would  require  the  NAT  =  1  line  of  input: 

10  1  (1.0, 0.0)  (0.0, 0.0)  0.2  . 


2.3  SUBROUTINE  WGEOM 


If  IRGM  =  0  and  INWR  =  1  (see  READ  1) ,  then 
and  attachment  geometry  will  be  specified  by  the 
WGEOM,  written  by  the  user,  and  which  has  the 
form : 


the  wire 
subroutine 
following 


PROGRAM  INPUT 


SUBROUTINE  WGEOM ( IA, IB , X , Y , Z , NM , NP , NAT, NSA, NPLA, VGA, 

2BDSK, ZLDA, NWG, VG, ZLD  WV) 

DIMENSION IA ( 1 ) ,IB(1),X(1),Y(1),Z(1) ,NSA(1) ,NPLA(1) ,BDSK(1) 
COMPLEXVGA(l) ,ZLDA(1) ,VG(1) ,ZLD(1) 


•  • 

MAIN  BODY  OF  SUBROUTINE 


RETURN 

END 


Before  describing  the  use  of  subroutine  WGEOM  in  more  detail 
we  will  simply  define  the  inputs  and  outputs  of  the  CALL 
statement . 

WGEOM  INPUTS: 

NWG  =  index  of  DO  faOO  loop  (see  Figure  1  and  READ  1). 

WV  =  wavelength  in  meters. 

Note  that  the  inputs  are  automatically  defined  in  the  MAIN 
program,  and  are  not  explicitly  defined  by  the  user.  They 
are  provided  as  an  aid  in  writing  WGEOM,  and  need  not  be 
used. 

WGEOM  OUTPUTS: 

I A  ( I )  =  endpoint  A  of  segment  I  (1  <=  I  <=  NM) . 

IB ( I )  =  endpoint  B  of  segment  I  (1  <=  I  <=  NM) . 

X ( J) ,Y (J) , Z ( J)  =  x,y,z  coordinates  in  meters  of  point 
J  (1  <=  J  <=  NP) . 

NM  =  total  number  of  wire  segments. 

NP  =  total  number  of  wire  points. 

NAT  =  total  number  of  attachment  points,  i.e.,  the  number 
of  points  at  which  wire  segments  contact  plates. 

NSA (K)  =  attachment  "location"  for  attachment  K 
(1  <=  K  <=  NAT) . 

NPLA (K)  =  plate  number  for  attachment  K  (1  <=  K  <=  NAT). 

VGA ( K )  =  complex  voltage  generator  (volts)  at  attachment 
K,  positive  polarity  points  up  wire  and  away  from 
plate  (1  <=  K  <=  NAT) 

BDSK(K)  =  disk  radius  for  attachment  K,  typically 
chosen  between  0.1  to  0.25  wavelength,  however, 
disk  must  not  overlap  an  edge  of  the  plate 
(1  <=  K  <=  NAT) . 

ZLDA (K)  =  complex  load  impedance  (ohms)  at  attachment 
K  (1  <=  K  <=  NAT) . 
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VG ( L)  =  complex  voltage  generator  (volts)  at  "location 
L"  in  the  wire  (1  <=  L  <=  2*NM) . 

ZLD(L)  =  complex  load  impedance  (ohms)  at  "location  L" 
in  the  wire  (1  <=  L  <=  2*NM) . 

All  of  the  above  outputs  must  be  defined  by  the  user  via 
FORTRAN  statements  in  subroutine  WGEOM . 

In  defining  the  arrays  NSA,  VG,  and  ZLD  we  referred  to 
a  "location"  in  the  wire.  This  means  either  by  point  A  or 
by  point  B  of  a  segment  of  the  wire  structure.  Specifically 
"location  L"  means: 

by  point  A  of  segment  L  if  L  <=  NM 

by  point  B  of  segment  NM-L  if  NM  <  L  <=  2*NM. 

For  example  ,  if  NM  =  8,  then  NSA(3)  =  1  means  that 
attachment  3  is  by  point  A  of  segment  1.  If  NM  =  8,  then 
NSA(3)  =  10  means  that  attachment  3  is  by  point  B  of  segment 
2.  If  NM  =  8,  then  VG(5)  =  (3. 0,2.0)  means  that  there  is  a 
generator  of  3.0  +  j 2 . 0  volts  by  point  A  of  segment  5, 
polarity  from  point  A  to  B.  If  NM  =  8,  then  ZLD(13)  = 
(1.0, -1.0)  means  that  there  is  a  load  of  1.0  -  jl.O  ohms  by 
point  B  of  segment  5. 

Note  that  only  the  non-zero  entries  in  VG  and  ZLD  need 
be  defined.  The  first  NAT  entries  in  NSA,  NPLA,  VGA,  BDSK , 
and  ZLDA  must  be  defined. 

After  writing  WGEOM,  the  user  has  two  choices  for 
including  it  as  part  of  the  code.  First  one  could  append  a 
FORTRAN  version  of  WGEOM  to  the  end  of  the  code.  This  has 
the  disadvantage  that  the  entire  code  must  be  compiled  each 
time  a  change  is  made  to  WGEOM.  A  preferred  method  would  to 
write  WGEOM  as  a  separate  file,  compile  it  to  form  an  object 
file,  and  link  the  object  file  with  an  object  file  of  the 
main  code. 

A  was  mentioned  above  it  is  advantageous  to  use 
subroutine  WGEOM  to  describe  the  wire  geometry,  rather  than 
the  input  file,  when  the  wire  structure  is  geometrically 
very  regular.  There  is  no  better  example  of  this  than  a 
straight  wire.  Thus,  suppose  one  wishes  to  study  the 
characteristics  of  a  center  fed  dipole  of  various  lengths 
and  segmentations.  Of  course  one  could  set  IRGM  =  1  in  READ 
1  and  input  each  new  geometry  via  READ'S  11-15.  The 
alternative  is  to  set  IRGM  =  0  and  write  a  subroutine  WGEOM 
capable  of  generating  the  dipole  geometry  for  arbitrary 
length  and  segmentation.  Figure  4  shows  a  straight  wire 
alligned  with  the  z-axis,  of  length  H,  and  divided  into  NM 
equal  segments.  We  note  the  following  concerning  this 
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00050 
00060 
00070 
00080 
00090 
00100 
00110 
00120 
00130 
00140 
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0017  0 
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00280 
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SUBKOUT1 NEWG EOM ( 1 A , I B , X , Y , Z ,  NM , NP, NAT, NSA , NPLA , VGA , BOSK , 

2  ZLPA , NWG , VG , ZLD , WV ) 

DIMENSIONIA(I) , 1B( 1 ) , X (1 ) , Y( 1 ) ,Z( 1 ) ,NSA(1 ) ,NPLA(1 ) , BDSK( 1 ) 
COMPLEXVGA(l) ,ZLDA(1) , VG( 1 ) ,7LD( 1 ) 

C 

C  GEOMETRY  FOR  A  CENTER  FED  DIPOLE. 

C 

C  SPECIFY  H  =  WIRE  LENGTH  AND  NM  =  NUMBER  OF  SEGMENTS. 

H=0 . 5 
NM-4 

C  INSURE  THAT  NM  IS  AN  EVEN  NUMBER. 

NM=2*( (NM+1) / 2) 

C  THE  NUMBER  OF  POINTS  IS 

NP=NM+1 

C  THE  SEGMENT  SIZE  IS 

DU=H/NM 

C  DEFINE  COORDINATES  OF  NP  POINTS  AND  NM  SEGMENTS. 

D01001=T , NP 
X(1)=0.0 
Y ( 1 ) =0  .0 
Z(  1  )  =  (  1-1 )  *DI1 
1F(  I .  F,Q.  NPIGOTOIOO 
IA  < 1 ) =1 
IB( 1 )=I+1 
100  CONTINUE 

C  DEFINE  GENERATOR  LOCATION  AND  VALUE. 

IGN-NM/2+1 

\G< IGN) -(1 .0,0.0) 

C  INDICATE  NO  ATTACHMENTS. 

NAT=0 

RETURN 

END 


Figure  5  -  A  subroutine  WGEOM  to  describe  the 
center  fed  dipole  of  Figure  4. 


PROGRAM  INPUT 


geometry : 

1.  The  number  of  points  is  NP  =  NM+1  . 

2.  The  segment  size  is  DH  =  H/NM  . 

3.  The  J-th  point  is  at: 

X(J)  =  0.0 
Y ( J )  =  0.0 
Z(J)  =  ( J-l) *DH  . 

4.  The  I-th  segment  has  endpoints: 

IA  ( I )  =  I 
IB { I )  =  1+1  . 

5.  If  the  dipole  is  to  be  center  fed  then  NM  must  be 
an  even  number,  and  the  generator  location  is: 

IGN  =  (NM/2) +1  or 

IGN  =  NM+NM/2  . 

Based  on  the  observations,  Figure  5  shows  a  subroutine  WGEOM 
which  generates  the  dipole  geometry.  COMMENT  statements 
describe  the  various  sections  of  the  subroutine.  As  shown 
in  Figure  5,  WGEOM  is  set  up  for  an  H  =  0.5  meter  dipole 
with  NM  =  4  segments.  The  advantage  of  writing  subroutine 
WGEOM  is  that  dipoles  of  different  lengths  and  segmentations 
can  be  obtained  by  simply  changing  lines  90  and  100.  Note 
that  excluding  COMMENT  statements,  the  main  body  of  this 
particular  WGEOM  contains  only  16  lines.  The  output  for 
this  geometry  will  be  shown  in  Design  Example  2  in  the  next 
chapter . 

As  a  second  example  of  writing  a  subroutine  WGEOM, 
consider  the  problem  of  describing  a  polygon  loop  of 
arbitrary  radius,  number  of  sides,  and  maximum  segment  size 
in  wavelengths.  We  will  use  the  notation: 

R  =  the  loop  radius  radius  in  meters. 

NS  =  the  number  of  sides  in  the  polygon  loop. 

SWX  =  the  maximum  segment  size  in  wavelengths. 

Figure  6  shows  a  hexagon  loop  with  NMS  =  2  segments  per 

side.  For  a  general  polygon  loop  we  note  the  following: 


PROGRAM  INPUT 


1.  If  there  are  NMS  segments  per  side,  then  the  total 
number  of  segments  is  NM  =  NMS*NS,  and  NP  =  NM. 

2.  If  SL  is  the  length  of  one  side  of  the  polygon, 
then  the  number  of  segments  per  side  is  the  first 
integer  >=  SL/(SWX*WV). 

3.  The  I-th  side  of  the  polygon  goes  from  phi  = 
( I-l ) *360/NS  to  I *360/NS  degrees. 

4.  Segment  J  has  endpoints: 

IA  ( J )  =  J 

IB ( J)  =  J+l  except  IB (NM)  =  1  . 

5.  If  the  loop  is  to  be  fed  at  phi  =  0.0,  then  the 
generator  location  is  at  IGN  =  1  or  IGN  =  2*NM. 


Using  these  observations,  Figure  7  shows  a  subroutine 
WGEOM  for  the  polygon  loop.  It  is  set  for  an  NS  =  6  sided 
loop  of  radius  R  =  0.3  meters  and  with  segments  less  than 
SWX  =  0.2  wavelengths  (regardless  of  the  frequency).  By 
changing  lines  100,  110,  and  120  one  can  easily  generate  the 
geometry  for  polygon  loops  of  different  radius,  number  of 
sides,  and  maximum  segment  size  in  wavelengths.  Note  that 
since  we  specified  the  maximum  segment  size  in  wavelengths, 
this  routine  can  be  used  at  virtually  any  frequency  with  no 
modification.  The  routine  automatically  increases  the 
number  of  segments  as  the  wavelength  decreases  (frequency 
increases)  and  decreases  the  number  of  segments  as  the 
wavelength  increases  (frequency  decreases) .  This  frequency 
independent  quality  is  especially  desireable  if  one  is  going 
to  analyze  the  wire  geometry  over  a  very  broad  frequency 
range. 
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00110 
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00140 
00150 
00160 
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00180 
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00270 
00280 
00290 
00300 
00310 
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00360 
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003  80 
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00400 
004  1  0 
00420 
00430 
00440 
00450 
00460 
00470 
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00490 
00500 
00510 
00520 
00530 
00540 
005  50 
00560 
00570 
00580 


SUBROUTI  NEW  EOM  ( I A ,  IB  ,1.1.1,  NM,  NP,  NAT,  NSA,  NPLA,  VGA,  BDSI, 

2  ZLDA ,  NWG ,  VG ,  ZLD ,  W  > 

D1  MENS  IONI  Ad  ) ,  IB(  1 ) , X ( I )  . Y ( 1)  ,2  ( 1 )  .NSA (I )  , NPL A d ) , BBS* d ) 
COMPLEX VG  A  ( I ) ,  7.LD A  ( 1  > .  VC  ( 1 ) ,  2  LD  ( I ) 

C 

C  GEOMETRY  FOR  POLYGON  LOOP. 

C 

C  SPECIFY  LOOP  R  -  LOOP  RADII’S  IN  METERS,  NS  =■  NUMBER  OF 

C  SIDES  IN  POLYGON  LOOP,  AND  SWI  -  MAXIMUM  SEGMENT  SIZE 

C  IN  RAVEL  ENG TBS. 

R-0 .3 
NS' 6 
SWX«0.2 

C  FIND  SL  *■  SIDE  LENGTH. 

PI  '4 ,0'ATANd  .0) 

DPII'2.0»P1/NS 

SL=2 ,0*R*SlN(DPH/2.0) 

C  FIND  NMS  -  NUMBER  OF  SEGMENTS  PER  SIDE  AND  DSL  '  THE 

C  SEGMENT  LENGTH . 

DSL' SVX'WV 
NMS'O. 99+51 /DSL 
DSL=SL/NMS 

C  FIND  NM  >  THF.  TOTAL  NUMBER  OF  SEGMENTS  AND  NP  -  THE  TOTAL 

C  NUMBER  OE  POINTS. 

NM'N\S«NMS 

NP=NM 

C  DEFINE  NMS  POINTS  AND  SEGMENTS  ON  EACH  OE  THE  NS  SIDES. 

DOIOOI'I ,NS 

C  THE  COORDINATES  OF  THE  FIRST  END  OF  SIDE  I  IS  AT 

PH1'( 1-1 ) *DPH 
Xl-R+COSU'Ul) 

Y1  =  R*S1N(PI11) 

C  THE  COORDINATES  OF  THE  SECOND  END  OF  SIDE  I  IS  AT 

PH2=1«DPU 
X2'R*COS(PU2) 

Y2»R*SIN<PH2) 

C  EACH  POINT  ON  SIDE  I  WILL  HE 

DX12'(X2-X1)/NMS 
DY12“(Y2-Y1)/NMS 

C  FROM  THE  LAST  POINT  ON  SIDE  I. 

DO200J-l,NMS 

C  DKFINF  THE  K  TH  POINT  AND  SEGMENT. 

R«(I-1)*NMS+J 

X(K)'X1+(J-I)#DX12 

Y(K)“Y1+(J-I)*DYI2 

Z(K)'0.0 

1A(K)'K 

1H(K  )  -K+  1 

1F(K.E0.NM)IB(C)'I 
200  CONTINUE 
100  CONTINUE 

C  PLACE  A  I  VOLT  GENERATOR  AT  THE  X  AXIS. 

IGN'I 

VG(IGN)'(I. 0,0.0) 

C  INDICATE  NO  ATTACHMENTS . 

NAT-0 

RETURN 

END 


Figure  7  -  A  subroutine  WGEOM  to  describe  the  polyqon  loop 
of  Figure  6. 


> 


OOIOO 

1211104  10  18  11 

00200 

1  1  3  0.0 

00300 

013  90.0 

00400 

0  1  3  0.0  90.0  0.0 

00500 

013  90.0 

00600 

150.0  38.0  0.001 

00700 

1  2 

00800 

4  4  3 

00900 

-0.5  -0.5  0.0 

01000 

0.5  -0.5  0.0 

01100 

0.5  0.5  0.0 

01200 

0  0 

01300 

3  4  11 

01400 

0.0  0.0  0.0 

01500 

0.0  0.0  0.25 

01600 

0.0  0.0  0.5 

01700 

-0.3  0.0  0.25 

01800 

1  2 

01900 

2  3 

02000 

2  4 

02100 

3  0  (0.0, 0.0)  (50.0,0.0) 

02200 

101  (1.0, 0.0)  (0.0, 0.0) 

Figure  8  -  Input  file  for  Example  1. 


CHAPTER  3 


DESIGN  EXAMPLES 


3.1  INTRODUCTION 

This  section  will  present  several  design  or  example 
runs  illustrating  the  use  of  the  code.  There  are  three 
purposes  to  thes  example  runs: 

1.  to  illustrate  input  data 

2.  to  illustrate  output  data 

3.  to  provide  trail  or  debugging  runs  for  a  new  user. 


3 . 2  EXAMPLE  1 

'*For  the  wire  and  plate  geometry  of  Figure  3  we  wish  to 
compute  the  currents,  impedance,  and  far-zone  elevation 
plane1  pattern  in  the  plane  phi  =  0.0.  The  wire  is 
considered  to  be  in  the  center  of  a  1  meter  square  plate. 
The  frequency  is  150  MHZ,  and  the  wire  is  taken  as  aluminum 
(conductivity  =  38  megamho/meter)  with  a  radius  of  0.001 
meter.  The  input  file  for  this  run  is  shown  in  Figure  8. 

A  three  view  orthagraphic  plot  of  the  wire  and  plate 
geometry  is  shown  in  Figure  9.  The  code  provides  this  plot 
if  NGO  =  0  in  READ  1.  Edges  of  the  plate  are  shown  as  solid 
lines.  Wire  segments  are  shown  as  solid  lines  with  small 
circles  at  the  endpoints.  A  summary  of  the  number  of  wire, 
plate,  and  attachment  modes  is  given,  as  well  as  a  scale 
indicating  what  one  inch  corresponds  to  in  wavelengths.  The 
plot  is  provided  for  two  reasons.  First  it  permits  the  user 
to  see  if  he  has  (probably)  specified  the  geometry 
correctly.  Secondly,  it  provides  a  permanent  pictoral 
documentation  of  the  geometry,  say  for  later  inclusion  in  a 
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report . 

The  output  for  this  run  is  shown  in  Appendix  I.  The 
file  is  called  OUTFL.DAT  on  the  ElectroScience  Lab  VAX 
system.  It  is  logical  unit  5.  The  output  begins  by  listing 
some  of  the  input  quanties  such  as  frequency,  wire  radius 
and  conductivity,  and  integration  parameters. 


Next  the  geometry  of  the  plates  is  specified.  Shown  is 
the  coordinates  of  three  consecutive  corners  of  each 
rectangular  plate,  together  with  the  plate  segmentation 
parameters  NM12,  NM23 ,  and  IPN.  If  IWR  =  1  in  READ  1,  then 
a  detailed  printout  of  the  overlapping  rectangular  surface 
patch  dipole  modes  is  provided.  Figure  10  shows  a  typical 
surface  patch  dipole  mode,  consisting  of  monopole  A  and 
monopole  B.  Monopole  A  is  defined  by  points  Al,  A2,  and  A3, 
while  monopole  B  is  define  by  points  B1 ,  B2,  and  B3 .  By 
convention  the  modal  current  flows  A2  to  Al  on  monopole  A, 
and  B1  to  B2  on  monopole  B.  If  IWR  =  1  (see  READ  1)  the 
coordinates  of  points  Al,  A2,  A3,  B1 ,  B2 ,  and  B3  are  printed 
for  each  surface  patch  mode. 

The  next  group  of  output  lists  the  wire  geometry. 

First  the  coordinates  in  meters  of  the  NP  points  are 
printed.  Next,  if  IWR  =  1,  the  wire  mode  layout  is  printed. 
Figure  11  shows  a  typical  wire  dipole  mode  going  from  point 
II  to  12  to  13  on  segments  JA  and  JB.  By  convention  the 
direction  of  positive  current  is  from  point  II  to  12  to  13. 

Next  the  endpoints  and  length  in  meters  of  each  of  the  NM 
segments  is  printed. 

The  next  group  of  output  is  th  attachment  geometry. 

For  each  attachment  the  following  is  listed: 

SEGMENT  =  wire  segment  which  contacts  the  plate. 

END  =  0  or  1  if  point  IA  or  IB  of  the  segment  contacts  the  plate. 
PLATE  =  plate  number  contacted  by  the  wire  segment. 

B  =  attachment  mode  disk  radius  in  meters. 

By  convention  the  attachment  mode  disk  current  flows  inward 
to  the  attachment  point  and  then  up  the  wire. 

The  next  printout  lists  the  complex  loads  and 
generators.  Following  this  NWR  =  the  number  of  wire  modes, 

NPLTM  =  the  number  of  plate  modes,  and  NAT  =  the  number  of 
attachment  modes  is  shown. 

The  output  terminates  at  this  point  if  NGO  =  0  (see 

READ  1) .  All  of  the  above  output  describes  the  detailed 
geometry  of  the  problem  specified  by  the  user.  The  user 
should  carefully  study  these  sections  to  make  sure  that  the 
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Figure  9  -  Three  view  plot  of  the  geometry  of  Example  1. 
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wire/plate  geometry  is  correct.  The  usual  procedure  for 
doing  this  is  as  follows.  When  a  new  geometry  is  being  set 
up,  the  initial  runs  are  made  with  NRUNS  =  0  (see  READ  1) . 
After  studying  the  printout  of  the  geometry  and  the 
orthographic  plot  of  the  geometry  and  being  convinced  that 

the  geometry  is  correct,  only  then  is  NGO  set  to  1  to 

perform  the  desired  computations. 

If  IWR  =  1,  the  next  set  of  output  lists  the  magnitude 
of  the  expansion  modes.  For  each  mode,  the  relative 
magnitude,  the  absolute  magnitude  (amps) ,  the  phase 
(degrees),  and  the  complex  magnitude  (amps)  is  tabulated. 
In  this  list  the  first  NWR  are  wire  modes,  the  next  NPLTM 
are  surface  patch  plate  modes,  and  the  last  NAT  modes  are 
attachment  modes. 

For  antenna  (as  opposed  to  scattering)  problems  the 
input  admittance  and  impedance  and  radiation  efficiency  are 
printed.  A  wire/plate  geometry  is  considered  to  be  a 

transmitting  antenna  if  ISE  =  0  in  READ  4  and  ISA  =  0  in 

READ  5.  THe  values  printed  for  input  admittance  and 
impedance  are  valid  if  there  is  only  one  generator  in  the 
wire,  and  it  is  1  +j0  volts. 

The  final  printout  is  the  far-zone  patterns.  For 
radiation  patterns,  the  gain  in  dB  for  theta  and  phi 
polarizations  is  shown.  For  scattering  patterns,  the  cross 
section,  sigma/wavelength**2  in  dB,  for  theta  and  phi  and 
cross  polarizations  is  shown.  The  phase  of  the  scattered 
field  field  is  also  shown  (the  incident  wave  has  zero  phase 
at  the  origon) .  For  backscatter ing  STPM  =  SPTM,  however, 
both  are  printed  as  a  check.  Due  to  lack  of  available  data 
for  comparison,  the  cross  polarized  data  is  not  as  reliable 
as  the  principle  polarization  data.  The  various  quanties 
are  defined  as  follows: 

STTM  =  scattering  cross  section  with  incident  and  scattered 
fields  theta  polarized. 

SPPM  =  scattering  cross  section  with  incident  and  scattered 
fields  phi  polarized. 

STPM  =  scattering  cross  section  with  incident  field  theta 
polarized  and  scattered  field  phi  polarized. 

SPTM  =  scattering  cross  section  with  incident  field  phi 
polarized  and  scattered  field  theta  polarized. 

GTHETA  =  gain  for  theta  polarization. 

GPHI  =  gain  for  phi  polarization. 


At  the  conclusion  of  each  run,  the  CPU  time  is  printed. 
In  this  case  is  about  117  seconds. 


I 


00100 

1  2 

1 

1  0  1  4  10  18  1 

00200 

0  1 

3 

0.0 

00300 

0  1 

3 

90.0 

00400 

0  1 

3 

0.0  90.0  0.0 

00500 

0  1 

3 

90.0 

00600 

300 

.0 

-1.0  0.001 

00700 

0  2 

00800 

0  0 

Figure  12  -  Input  file  for  Example  2. 
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3.3  EXAMPLE  2 

Here  we  wish  to  find  the  input  impedance  of  a  center 
fed  half-wave  dipole  of  radius  0.001  wavelength,  and 
assuming  perfect  conductivity  of  the  wire.  The  subroutine 
WGEOM,  shown  in  Figure  5  sets  up  a  center  fed  dipole  of 
length  0.5  meters  and  with  4  segments.  It  will  be  a 
half-wave  dipole  if  the  frequency  =  300  MHZ.  Figure  12 
shows  the  required  input  file  for  this  problem.  Note  that 
in  READ  1,  IWRZT  =  1  which  will  result  in  a  printout  of  the 
moment  method  impedance  matrix.  Also,  in  READ  1  IRGM  =  0, 
so  that  the  wire  geometry  is  to  be  generated  in  subroutine 
WGEOM  (rather  that  from  READ'S  11-15  of  the  input  file).  A 
three  view  orthographic  plot  of  the  dipole  is  shown  in 
Figure  13. 

The  output  is  shown  in  Appendix  II.  The  only  portion 
of  the  output  not  previously  described  is  the  printing  of 
the  impedance  matrix.  This  occurs  just  above  the  printout 
of  the  input  admittance.  The  lower  triangular  part  of  the 
symmetric  impedance  matrix  is  printed  by  columns.  Following 
this,  the  input  impedance  of  the  half-wave  dipole  is  shown 
as  about  81  +j41  ohms.  Note  that  this  run  took  about  0.5 
seconds . 


3.4  EXAMPLE  3 

In  this  example,  we  will  illustrate  a  scattering 
computation  involving  intersecting  plates.  Say,  for 
example,  that  we  wish  to  compute  the  backscatter  from  the 
corner  reflector  shown  in  Figure  14.  It  consists  of  two  1.0 
by  0.5  wavelength  plates,  intersecting  along  the  z  axis. 
The  input  for  this  geometry  is  shown  in  Figure  15.  A 
backscatter  azimuth  pattern  in  the  plane  theta  =  90.0 
degrees  is  specified  by  setting  ISA  =  1  and  THSA  =  90.0  in 
READ  5.  In  READ  7  IOVL  =  2  so  that  overlap  modes  will  be 
inserted  to  connect  the  two  plates.  Figure  16  show  the 
orthagraphic  view  and  Appendix  III  shows  the  output.  Note 
that  after  specifying  the  coordinates  of  plates  1  and  2,  the 
output  indicates  that  4  overlap  modes  were  inserted 
connecting  side  1  of  plate  1  to  side  1  of  plate  2.  See 
Figure  2  for  definition  of  plate  sides.  The  backscatter 
patterns  for  theta  (STTM)  and  phi  (SPPM)  polarization  are 
shown,  and  plotted  in  Figures  17a, b.  The  cross  polarization 
is  not  plotted  since  it  is  so  small. 
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7a  -  Theta  polarized  azimuth  backscatter 
pattern  for  Example  3. 
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3.5  EXAMPLE  4 


Example  4  is  identical  to  Example  3  except  that  we  wish 
a  bistatic  scattering  pattern,  with  the  incident  wave 
comming  from  theta  =  90.0  degrees  and  phi  =  45.0  degrees. 
Thus,  we  now  set  THIN  =  90.0  and  PHIN  =  45.0  in  READ  4,  and 
ISA  =  2  in  READ  5.  The  input  file  is  shown  in  Figure  18, 
and  Appendix  IV  shows  the  output.  The  bistatic  scattering 
patterns  are  shown  in  Figures  19a, b. 


3.6  EXAMPLE  5 

This  example  will  illustrate  the  use  of  READ  10  to  save 
computational  time.  Say  we  are  given  the  problem  of 
computing  the  input  impedance  of  a  quarter  wave  monopole  at 
several  locations  on  plate  2  of  the  three  plate  bend 
illustrated  in  Figure  20a.  Here  we  will  compute  the  input 
impedance  for  the  monopole  attached  at  location  1  =  (x,y,z) 
=  (0 .0 ,0 .0 ,0 .0)  and  at  location  2  =  (x,y,z)  =  ( 0 . 0 , 0 . 3 , 0 . 0 )  . 
Computation  time  will  be  saved  by  noting  that  when  the 
monopole  moves  from  location  1  to  location  2,  the  plate 
geometry  does  not  change.  Thus  all  of  the  plate  to  plate 
impedances  do  not  change. 

To  find  the  impedance  at  these  two  locations  we  will 
set  NWG  =  2  in  READ  1,  indicating  that  we  are  running  two 
tfire  geometries.  For  the  first  geometry  we  will  set  IWRZM  = 
1  and  IRrZM  =  0  in  READ  10.  Thus  for  the  first  geometry  we 
will  compute  the  entire  impedance  matrix  and  then  write  it, 
unformatted,  onto  the  disk  file  ZMAT.DAT. 1  (logical  unit  1). 
For  the  second  geometry  we  will  set  IRDZM  =  2,  indicating 
that  the  impedance  matrix  is  to  be  read  in,  unformatted, 
from  file  ZMAT. DAT. 1  on  logical  unit  1.  Also  the  plate  to 
plate  block  of  the  impedance  matrix  is  not  be  recomputed, 
resulting  in  a  savings  of  time.  For  the  second  geometry 
IWRZT  may  be  set  to  0  or  1.  Figure  21  shows  the  input  file. 
READ  10  is  at  lines  2000  and  2800.  The  monopole  is  composed 
of  two  segments  as  seen  in  Figure  20b.  The  wire  is 
perfectly  conducting  with  radius  of  0.001  meters.  Figure  22 
shows  the  three  view  orthographic  plot  with  the  monopole  at 
location  1.  the  output  file  is  shown  in  Appendix  V.  Note 
that  for  the  first  geometry  where  the  entire  impedance 
matrix  was  computed  (IRDZM  =  0)the  run  time  was  about  411 
seconds.  However,  for  the  second  geometry,  where  where  the 
impedance  matrix  was  read  in  and  the  plate  to  plate  block 
was  not  computed  (IRDZM  =  2),  the  time  was  reduced  to  98 
seconds. 


00100  1211004  10  18  01 

00200  0  1  3  0.0 

00300  0  1  3  90.0 

00400  0  1  3  0.0  90.0  45.0 

00500  2  1  5  90.0 

00600  300.0  -1.0  0.001 

00700  2  2 

00800  423 

00900  0.0  0.0  0.0 

01000  0.0  0.0  1.0 

01100  0.5  0.0  1.0 

01200  423 

01300  0.0  0.0  0.0 

01400  0.0  0.0  1.0 

01500  0.0  0.5  1.0 

01600  0  0 


Figure  18  -  Input  file  for  Example  4. 
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Figure  20a  -  Geometry  for  the  three  plate  hend 
of  Example  5.(b)  Geometry  for 
wire  of  Example  5 


00100  1212004  10  18  11 

00200  0  0  3  0.0 

00300  0  1  3  90.0 

00400  0  1  3  0.0  90.0  0.0 

00500  013  90.0 

00600  300.0  -1.0  0.001 

00700  3  2 

00800  443 

00900  0.5  -0.5  -1.0 

01000  0.5  -0.5  0.0 

01100  -0.5  -0.5  0.0 

01200  443 

01300  0.5  -0.5  0.0 

01400  0.5  0.5  0.0 

01500  -0.5  0.5  0.0 

01600  243 

01700  0.25  0.5  0.0 

01800  -0.25  0.5  0.0 

01900  -0.25  0.5  1.0 

02000  1  0 

02100  2310 

02200  0.0  0.0  0.0 

02300  0.0  0.0  0.125 

02400  0.0  0.0  0.25 

02500  1  2 

02600  2  3 

02700  1  0  2  (1.0, 0.0)  (0.0, 0.0)  0.2 

02800  0  2 

02900  2310 

03000  0.0  0.3  0.0 

03100  0.0  0.3  0.125 

03200  0.0  0.3  0.25 

03300  1  2 

03400  2  3 

03500  102  (1.0, 0.0)  (0.0, 0.0)  0.2 


Figure  21  -  Input  file  for  Example  5. 
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Figure  22  -  A  three  view  plot  of  the  geometry 
of  Example  5. 


CHAPTER  4 


ARRAY  DIMENSIONS  AND  FILE  DESCRIPTIONS 


4.1  ARRAY  DIMENSIONS 

The  array  dimensions  are  defined  by  DIMENSION  and 
COMPLEX  statements  located  near  the  top  of  the  main  program. 
All  arrays  have  either  fixed  dimensions,  independent  of  the 
geometry  being  run,  or  are  dimensioned  according  to  one  of 
the  following  dimension  indicators: 

INM  =  maximum  number  of  wire  segments. 

ICJ  =  maximum  number  of  wire  modes. 

IDC  =  maximum  number  of  elements  in  wire/wire  block  of 

impedance  matrix. 

IPLM  =  maximum  number  of  plate  modes. 

IPL  =  maximum  number  of  plates. 

I AT  =  maximum  number  number  of  wire  to  plate  attachments. 

INP  =  maximum  number  of  wire  points. 

ITOT  =  maximum  total  number  of  modes  (wire  +  plate  +  attach.) 
IDZT  =  maximum  number  of  elements  in  impedance  matrix. 

The  dimension  indicators  are  defined  below  the  DIMENSION  and 
COMPLEX  statements,  and  typically  have  the  values: 

INM  =  491 
ICJ  =  492 

IDC  =  ( ICJ* ICJ + ICJ ) / 2 

IPLM  =  490 

IPL  =  7 

I  AT  =  2 

INP  =  493 

ITOT  =495 

IDZT  =  ( ITOT*ITOT+ITOT) / 2 

Thus  the  program  is  typically  set  for  no  more  than  INM  =  491 
wire  segments,  ICJ  =  492  wire  modes,  IDC  =  121,278  elements 
in  the  (lower  triangular  part  of  the)  wire  to  wire  block  of 
the  impedance  matrix,  IPLM  =  490  plate  modes,  IPL  =  7 
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plates,  I AT  =  2  attachment  points  or  modes,  INP  =  493  wire 

points,  ITOT  =  495  total  modes,  and  IDZT  =  122,760  elements 
in  the  (lower  triangular  part  of  the  symmetric)  impedance 
matrix.  Note  that  while  the  number  wire  modes  can  be  up  to 
492,  the  number  of  plate  modes  up  to  490,  or  the  number  of 
attachment  modes  up  to  2,  the  total  number  of  modes  can  not 
exceed  495. 

There  are  two  steps  required  to  change  the  allowed 
dimensions : 

1.  One  must  change  the  appropiate  dimension  indicator. 

2.  One  must  re-dimension  all  arrays  associated  with 
that  dimension  indicator.  Arrays  dimensioned  by 
the  same  indicator  are  grouped  together,  and  are 
clearly  marked  by  COMMENT  statements. 


4.2  FILE  DESCRIPTIONS 

A  description  of  program  or  source  files  and  I/O  files 
will  now  be  given. 


4.2.1  PROGRAM  FILES 

The  computer  code  is  contained  in  several  different 
(disk)  files  on  the  ElectroScience  Lab's  Digital  Equipment 
Corporation  VAX  computer  system.  A  listing  of  these  files 
follows.  The  file  types  are  shown  as  FORTRAN  (except  for 
'PLOTLIB  which  an  object  file). 

STDMM . FOR  -  the  main  program  plus  various  surface  patch 
subroutines 

THNWRS.FOR  -  various  thin  wire  subroutines 

WGEOM. FOR  -  subroutine,  written  by  the  user,  describing  the 
WIRE  geometry  (see  section  2.3) 

LIB. FOR  -  various  special  library  routines.  At  present  we 
use  only  the  function  subroutine  GETCP(I)  where  I  =  clock 
reading  in  hundreths  of  a  second.  Since  this  routine  tends 
to  be  hardware  dependent,  it  is  not  included  when  the  program 
is  sent  outside  the  ESL. 

GPLOT. FOR  -  subroutine  to  a  make  three  view  orthographic  plot  of 
wire  and  plate  geometry 

'PLOTLIB  -  contains  various  plotting  subroutines.  If  the  code  is 
being  sent  outside  the  ESL  some  routines  must  be  omitted  due  to 
contactural  restrictions.  When  this  file  is  supplied  to  an 
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outside  user  it  will  be  termed  PLOTLIB.FOR.  Of  the  many 
subroutines  omitted, the  only  four  used  in  the  code  are: 

VPLOTS (1,0 ,0) 

VPLOTS  reserves  the  plotter. 

1=1  implies  the  plot  is  for  the  Versatek  paper  plotter 
2  implies  the  plot  is  for  the  Megatek  CRT  plotter 
0  implies  the  program  gives  the  user  a  choice  of  plotter 

PLOT (X , Y , I ) 

PLOT  moves  the  the  pen,  with  pen  up  or  down. 

X,Y  implies  move  pen  to  these  coordinates  (inches) 

1=2  implies  lower  pen  before  moving 

=  3  implies  raise  pen  before  moving 

=  -2  or  -3  implies  same  as  2  or  3  except  reset  origon  after 
moving 

=  -999  implies  go  to  lower  left  corner  of  next  page  with  pen 
up  and  reset  origon 

=  999  implies  this  is  the  last  plotting  call  so  plot  everything 
in  the  plot  buffer  and  release  the  plotter 

A  call  to  VPLOTS  must  precede  all  plotting,  and  call  to  PLOT  with 

I  =  999  must  be  the  last  plotting  call. 

NU MB ER(X,Y,HT,FPN, ANGLE, N) 

NUMBER  plots  out  a  floating  point  number. 

X,Y  =  coordinates  of  lower  corner  of  output  number  in  inches 

HT  =  height  of  the  output  number  in  inches.  If  HT  >  0,  then  the 

output  will  be  plotted  to  the  right  of  X,Y?  if  HT  <  0,  it  will 

be  plotted  to  the  left. 

FPN  =  floating  point  number  to  be  plotted 

ANGLE  =  angle  in  degrees  (counterclockwise)  with  respect  to  the 
X  axis  at  which  the  output  number  is  to  be  plotted 

N  =  an  integer  specifying  the  output  format.  If  ABS(N)  <  100,  then 

FPN  will  be  plotted  in  the  "F"  format.  If  N  >  0,  then  N  digits 

will  be  plotted  after  the  decimal  point  in  addition  to  all  digits 
before  the  decimal  point.  If  N  <  0,  then  no  digits  will  be  plotted 
after  the  decimal  point,  and  the  decimal  point  plus  the  first 
-N-l  digits  to  the  left  of  the  decimal  point  will  be  suppressed. 

If  ABS(N)  >=  100,  then  FPN  will  be  plotted  in  an  E  format, 
that  is,  the  mantissa  of  FPN  followed  by  an  "X",  followed  by  a  10 
raised  to  a  power.  If  N  >  100,  then  there  will  be  one  digit  to  the 
left  and  N-100  digits  to  the  right  of  the  decimal  point  in  the 
mantissa.  If  N  <  -100,  then  an  integer  mantissa  of  -N-100  digits 
followed  by  a  power  of  10  will  be  plotted. 

SYMBOL ( X , Y , HT , LAB  EL , ANGL  E , NC ) 

SYMBOL  plots  a  character  or  string  of  characters. 
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X,Y  =  coordinates  in  inches  of  the  lower  left  hand  corner  of  the 
symbol  to  be  drawn 

HT  =  the  height  in  inches  of  the  character  to  be  drawn.  HT  should 
be  a  multiple  of  7  times  the  plotter  increment. 

LABEL  =  if  NC  >  0,  then  LABEL  is  a  literal  variable  or  constant 

representing  the  character  string  to  be  plotted.  NC  =  the  number 
of  characters  to  be  plotted. 

=  if  NC  =  -1,  then  LABEL  is  an  integer  expression  ranging  from 
0  to  127  which  represents  a  single  character  to  be  plotted.  These 
symbols  and  their  codes  are  given  in  Figure  23. 

ANGLE  =  angle  in  degrees  between  the  symbol  to  be  plotted  and  the 
X  axis 

NC  =  see  LABEL 


If  the  user  can  not  supply  a  subroutine  GETCP,  then  all 
references  to  this  subroutine  must  be  deleted,  and  the 
program  will  not  supply  run  time  information.  If 
subroutines  NUMBER  SYMBOL  VPLOTS  and  PLOT  are  not  available, 
then  all  calls  to  plotting  subroutines  must  be  eliminated. 
It  is  felt  in  most  cases,  if  a  plotter  is  available,  then 
these  routines  will  also  be  available.  If  plotting  is  not 
desired  then  subroutines  POLAR  and  GPLOT  should  be  removed. 
All  calls  to  these  and  to  subroutines  PLOT,  VLOTS,  and 
SYMBOL  should  be  removed  or  COMMENTED  out. 

In  summary  when  the  code  is  supplied  outside  the  ESL 
the  following  FORTRAN  files  are  included  in  a  single  file 
called  OSUESP . FOR: 

STDMM. FOR 
THNWRS. FOR 

WGEOM. FOR  (see  Figure  5) 

GPLOT. FOR 
PLOTLIB . FOR 

I-  te  that  the  subroutine  WGEOM  supplied  is  for  a  dipole.  To 
obtain  a  new  geometry,  then  he  must  write  a  new  subroutine 
WGEOM  and  replace  the  one  supplied  (see  Chapter  2  on  READ  1 
and  SUBROUTINE  WGEOM) .  Unless  otherwise  indicated  the 
magnetic  tape  format  is: 


ANSI  ASCII  label  -  OSUESP 
9  track  800  bits/inch 
Not  blocked  (80  characters/record) 
ASCII  character  set 
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Figure  23  -  Symbol  table. 
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4.2.2  INPUT/OUTPUT  FILES 

The  following  files  are  used  to  input  and  output  data. 

INSTD.DAT  -  the  input  file  described  in  Chapter  2.  It  must  be 
defined  as  logical  unit  5  by  the  user. 

OUTFL.DAT  -  the  output  file  as  described  in  Chapter  3  and  shown 
in  the  appendices.  It  is  defined  as  logical  unit  6  via 
OPEN  and  CLOSE  statement  in  the  main  program. 

ZMAT.DAT  -  input  and  output  of  the  impedance  matrix  as  described 
in  section  2.2.5.  It  is  defined  as  logical  unit  1  via  OPEN 
and  CLOSE  statements  in  the  main  program. 
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o 


a.  r<#NOK1«CPPJtf10Cr<»Pr'*r'v£P(\lf  (r»«*r-C,>^^P('jirff'-*PtJlrl>J)PrJira»«*NCj>Ov|  P  r  W  (T—  «P-0»0 
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I  JldKtX3  NOIS>0  Mi  indLClO 


OUTPUT  FOP  OF  SIGN  EXAMPLE  2  P*qe 


OUTPUT  FO*  OESIGN  C  X  4UPLF 


LISTING  OF  LOAOS  AND  GENFRATORS 


Cl  C*  *0  »«  |  T  #  U  iT  ^«  •»  *>•  vC  i,T  1/  if  i,"  ^  lT  •o  fv 

i  i  i  i  i  i—  ——  —  —  —  —  —  —  ~  —  ~  —  —  —  —  —  ~  — 

ii  l  i  i  i  i  i  i  i  i  i  ill 


CL  a  ^  N  ^  Lf>  Cl  c  ^  c.  ft  . 

in  ^  \T  lT  «r  ♦  k-  k-  K-  r  . 

I  I  I  I  I  I  I  I  I  I  I  I  I 
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—  —  I  I  —  ~  —  —  llllll 
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-*r'fvr»*|/*lT 


>r  o  «f  ir  —  \f.  ^  —  (N(»'o<r>^(T^cv--»ir'-- 

r-,O'**in#(,>fr'DrvrN'''fiPr'*r''<^0  p  •.  m  »  n  r  r  o  r-  r-  *■  o  n  *» 

y*  *  r  —  e>  c  \o^^^v-*'*cp«-'C,‘»ir,^occpM  —  c'C'<r«-of*«»c.  p  if  n  o  s 
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C'C'OOOOOoOl 
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0U1PUI  F3R  OCSIGN  £<A*PlC  3  P j-je 
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f^fr^cjirarfrirr^^irocrir  irin^iotrvrinror  <r  r  tro-T-oor-o^a 
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H^#triro^'rHNs»<rtfOK  f-  ff>  it  -  N(ro#ff'vfoi'*(rM»vfCT,  tfoc« 
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i  I  i  i  l  i 


rn»f\fNP'i,"ctp  iTC'fftr'ttfffjMrNCiri*'®  r,  er*-*-**e>er  o  r-  c  vf  r  ,  o 
oor  a  a  «r  «■  ff'vT'^nnovrrr.  crcc  —  rj^-vr^^c^cca:^  — 

<\i  ►-  f-  #•  r  r\fNir^(/lKvro'rj#s»‘Mr^(rpr  o  sr  «\i  it  m  —  it  r,  r-  r  ^  o 


'  ♦  ^  ^  «  O  H  ( 


■  n  n  h  c  -  f-  k>  #  r  v*i  ip  ♦  #  • 


■  ►«  «?  —  —  rj  ♦  it 


cfi^^ir-oHcrevf  «  ^ki»  ♦  ir»eoo*-ocie>*  o  *a  **  *»  <r  w  +  <r  <r  ~  sc 

ooravT^^HvTf  vfHNMTiwfNjrG^r'ti^s^S^CCf'^Cor^r  — 

tr  rj  ♦  ♦  r  ocjC'f'vtrf\inr\,*-<j-<\/N£cr<rv£r'j»?'C,-f'  fso^irooHrKio'o 
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0000000000004 
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>000000000000000000000000000000000000 


OUTPUT  FOR  DESIGN  f X AMPLE 


nur  =  numoer  of  u i re  moses  - 

A) PL  T  M  =  NIWBER  OF  PLATE  MODES  r 
NAT  =  NUMBER  0=  ATTACHMENT  MOOES 


Miijo  wti  indino 


■  — *  iT  tr  r 

-  P-  vT  Jl  C3  « 

I  I  I  —  I 


vT  —  iT  o  -'  p*  «»  tr  —  >->  f.  v.~  ^  f*  it  #  r  r“'  ^  ♦  ri  ^  T1  *■  iT  >'  n  <**-►-  l"  r  tr  ^  ^ 

«■  ^  f:  <\'  —  o  o  o  p-  *'"'  i  m  r  "  •  kfi  »p  «;  *a  *r  ip  it  # 

~  i  l  i  i  i  I  i  i  i  i  l  i  l  i  i  i  i  i  i  i  l  i  i  i  i  i  i 
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(\o^»ro^^*acQrr'CT''*,«noocr(^'*a'or^~*“\/itr^,^'fv»-*^  i  t\;  *^  *'■'  oi  i  «-  cm  *p 

•  •  *•••  •  •■•  i«  •  i  •  i  •  i  •  •  i  •  i  i  i  ( 

i  i  i 


k*  cc  (rrj°,'oririrj"'Mf  r-  of  c~  rr  p  «r  »  ir*  r-  ►-  tr-  f  o  ♦  ^  f  b  •  r  r  it  rvi  ?• 

c  H(t»ff'rt">trirr-''**KO(rt»‘OH(>.  *r  »*•  r  -.ir(^^-^rjf''r,,-'-'>’rrrfi'Xiirr 
•^ir-'c'iMtM^rrr  o  »r  i-  *•  f  f  cj  f  ♦  t"  ,*  o  »p  o 


*.r-.*-.p-.\fp-p-p- 

i  i  i  i  i  i  i  i 


<j,(ro#rHirf'irrsjin'>D‘-*'KX(r(rirNslf!r»''t^^''# 

Nsr(r-r(rsKN>.^^KKKKKKSSSKS».*.k.K*.KK 


h  it  f  »  p,  <»  —  tr  » 

r^Pv.rjtr^'fJiro  —  —  tr— ••-vc^-oir  <\  * 

<va:m^*«-orgo^cD^ip»'  k  ^  >cvP  -  —  « 


-P  «•  vf  ^ 

f-a  tfp.pciroi^poffMtfc  o  —  ♦  tr  «r 
otro<MO«'*OL"iKoc*»iro«rc\((N)t/'  r  it  <r 


p»#*c#ii'v:s(r  -  fSirMPN»fp*»»nir#*irtfh(co-»d  cctrhcifff#* 
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cr  Hr-i'N?‘N.1;(T>'i/'i"rwpff  <\CStfpNKfira  if  ♦  c  it  «■  <*“  ^  in  r  T  l'  r.  « 
—  tr  r  »r>ff‘fs.^-\fc:o'^p«-a:vrK'K.tr\f-*fv.cP.pjyc’  —  --►-♦♦vr»r«nsrtrr3a:ir«-'^ 
♦  c^ircoKur^h.trtprcMrr'tr.  fvfir^erL*-'a:*o^^r(''r'»'Hir#ov.',r  tr 


>  o  ~  r  p.  r,  ft  r.  r  r  rv.  r  pi  cm  pj  cm  <m  r.  pj  p.  —  —  . 


♦  »fl  »  P  PJ  P  r  P  CC  N  H  vA""^  f\DlPPPlP\COlPff^PPvfi^\flO|Pi.nirP(C(rP.lP»d*ce 

p.o«'p.P.«-«^-\rvrotr^-p  ^-cr  cp  it  \r  --  **«  *■  p  -  r  it  -  o  p  ^  vr  if  it  □  p 

##iPwOpH»dPvfPincc8:riPHvrp-»Hff'\fr<iP#pPro\1ipH^'0\f<rf'i'o# 
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01P  oiroiroiToiroiT  oiromoiroiroirorpirpif  o  -n  o  iP  o  o  i  -  o  it  o  tr  o 
iCvflsstrirr  r-o^-~-t\r„^»o<r^irw'vr«''p--p-a:!3‘OOr-io~-~-r  tM  ^  ♦ir*‘,vr 

^  H  H  H  H  K  H  H  PI  pi  Pi  P'  P.  n  (M  P.  PgPgfV(VP  ^Pjf;CVPy|p, 
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OUTPUT  FOR  DESIGN  EXAMPLE 


PLATE  NUMBER  T 
NMI2  -  2  NM2J  -  H  IP  -  1 

X,Y*Z  COOR.C METERS >  OF  CORNER  1  =  0.250  0.530  0.000 
X  *  Y  *  Z  COOR. (METERS)  OF  CORNER  2  =  -0.250  0.500  0.000 
X  *  Y  *  Z  C  OOR • ( ME  TrR3  >  OF  CORNER  3  =  -0.250  0.500  1.000 


GfOMETRY  FOR  THE  1  ATTACHMENT  POINT 


LISTING  OP  LOAOS  AND  GENERATORS 


PLATT  NUNRTR  ? 

\“1?  :  A  N-?  3  =  A  IP  =  3 

I,T*7  CCO^.CMETERSI  OF  COPNF*  1  r  ’'.SOO  -0.511  0.310 

COIR.  C**:fEa->  OF  CORNER  **  '  3.500  0.501  0.000 

A  •  T  *  ?  eqoR. C METERS )  OF  CORNER  3  =  -0.500  0.503  0.000 


COO*. t ME TE OF  24  MOOES  ON  THIS  PLATE 


;ESMrMS  ON  THE  WIPE 

u(j>  n(j)  o<jm«i 


OUTPUT  FOI  DESIGN  EXAMPLE  5  Page 


